3.5.57 \(\int \frac {1}{c+\frac {a}{x^6}+\frac {b}{x^3}} \, dx\) [457]

Optimal. Leaf size=631 \[ \frac {x}{c}+\frac {\left (b-\frac {b^2-2 a c}{\sqrt {b^2-4 a c}}\right ) \tan ^{-1}\left (\frac {1-\frac {2 \sqrt [3]{2} \sqrt [3]{c} x}{\sqrt [3]{b-\sqrt {b^2-4 a c}}}}{\sqrt {3}}\right )}{\sqrt [3]{2} \sqrt {3} c^{4/3} \left (b-\sqrt {b^2-4 a c}\right )^{2/3}}+\frac {\left (b+\frac {b^2-2 a c}{\sqrt {b^2-4 a c}}\right ) \tan ^{-1}\left (\frac {1-\frac {2 \sqrt [3]{2} \sqrt [3]{c} x}{\sqrt [3]{b+\sqrt {b^2-4 a c}}}}{\sqrt {3}}\right )}{\sqrt [3]{2} \sqrt {3} c^{4/3} \left (b+\sqrt {b^2-4 a c}\right )^{2/3}}-\frac {\left (b-\frac {b^2-2 a c}{\sqrt {b^2-4 a c}}\right ) \log \left (\sqrt [3]{b-\sqrt {b^2-4 a c}}+\sqrt [3]{2} \sqrt [3]{c} x\right )}{3 \sqrt [3]{2} c^{4/3} \left (b-\sqrt {b^2-4 a c}\right )^{2/3}}-\frac {\left (b+\frac {b^2-2 a c}{\sqrt {b^2-4 a c}}\right ) \log \left (\sqrt [3]{b+\sqrt {b^2-4 a c}}+\sqrt [3]{2} \sqrt [3]{c} x\right )}{3 \sqrt [3]{2} c^{4/3} \left (b+\sqrt {b^2-4 a c}\right )^{2/3}}+\frac {\left (b-\frac {b^2-2 a c}{\sqrt {b^2-4 a c}}\right ) \log \left (\left (b-\sqrt {b^2-4 a c}\right )^{2/3}-\sqrt [3]{2} \sqrt [3]{c} \sqrt [3]{b-\sqrt {b^2-4 a c}} x+2^{2/3} c^{2/3} x^2\right )}{6 \sqrt [3]{2} c^{4/3} \left (b-\sqrt {b^2-4 a c}\right )^{2/3}}+\frac {\left (b+\frac {b^2-2 a c}{\sqrt {b^2-4 a c}}\right ) \log \left (\left (b+\sqrt {b^2-4 a c}\right )^{2/3}-\sqrt [3]{2} \sqrt [3]{c} \sqrt [3]{b+\sqrt {b^2-4 a c}} x+2^{2/3} c^{2/3} x^2\right )}{6 \sqrt [3]{2} c^{4/3} \left (b+\sqrt {b^2-4 a c}\right )^{2/3}} \]

[Out]

x/c-1/6*ln(2^(1/3)*c^(1/3)*x+(b-(-4*a*c+b^2)^(1/2))^(1/3))*(b+(2*a*c-b^2)/(-4*a*c+b^2)^(1/2))*2^(2/3)/c^(4/3)/
(b-(-4*a*c+b^2)^(1/2))^(2/3)+1/12*ln(2^(2/3)*c^(2/3)*x^2-2^(1/3)*c^(1/3)*x*(b-(-4*a*c+b^2)^(1/2))^(1/3)+(b-(-4
*a*c+b^2)^(1/2))^(2/3))*(b+(2*a*c-b^2)/(-4*a*c+b^2)^(1/2))*2^(2/3)/c^(4/3)/(b-(-4*a*c+b^2)^(1/2))^(2/3)+1/6*ar
ctan(1/3*(1-2*2^(1/3)*c^(1/3)*x/(b-(-4*a*c+b^2)^(1/2))^(1/3))*3^(1/2))*(b+(2*a*c-b^2)/(-4*a*c+b^2)^(1/2))*2^(2
/3)/c^(4/3)*3^(1/2)/(b-(-4*a*c+b^2)^(1/2))^(2/3)-1/6*ln(2^(1/3)*c^(1/3)*x+(b+(-4*a*c+b^2)^(1/2))^(1/3))*(b+(-2
*a*c+b^2)/(-4*a*c+b^2)^(1/2))*2^(2/3)/c^(4/3)/(b+(-4*a*c+b^2)^(1/2))^(2/3)+1/12*ln(2^(2/3)*c^(2/3)*x^2-2^(1/3)
*c^(1/3)*x*(b+(-4*a*c+b^2)^(1/2))^(1/3)+(b+(-4*a*c+b^2)^(1/2))^(2/3))*(b+(-2*a*c+b^2)/(-4*a*c+b^2)^(1/2))*2^(2
/3)/c^(4/3)/(b+(-4*a*c+b^2)^(1/2))^(2/3)+1/6*arctan(1/3*(1-2*2^(1/3)*c^(1/3)*x/(b+(-4*a*c+b^2)^(1/2))^(1/3))*3
^(1/2))*(b+(-2*a*c+b^2)/(-4*a*c+b^2)^(1/2))*2^(2/3)/c^(4/3)*3^(1/2)/(b+(-4*a*c+b^2)^(1/2))^(2/3)

________________________________________________________________________________________

Rubi [A]
time = 0.79, antiderivative size = 631, normalized size of antiderivative = 1.00, number of steps used = 15, number of rules used = 9, integrand size = 14, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.643, Rules used = {1354, 1381, 1436, 206, 31, 648, 631, 210, 642} \begin {gather*} \frac {\left (b-\frac {b^2-2 a c}{\sqrt {b^2-4 a c}}\right ) \text {ArcTan}\left (\frac {1-\frac {2 \sqrt [3]{2} \sqrt [3]{c} x}{\sqrt [3]{b-\sqrt {b^2-4 a c}}}}{\sqrt {3}}\right )}{\sqrt [3]{2} \sqrt {3} c^{4/3} \left (b-\sqrt {b^2-4 a c}\right )^{2/3}}+\frac {\left (\frac {b^2-2 a c}{\sqrt {b^2-4 a c}}+b\right ) \text {ArcTan}\left (\frac {1-\frac {2 \sqrt [3]{2} \sqrt [3]{c} x}{\sqrt [3]{\sqrt {b^2-4 a c}+b}}}{\sqrt {3}}\right )}{\sqrt [3]{2} \sqrt {3} c^{4/3} \left (\sqrt {b^2-4 a c}+b\right )^{2/3}}+\frac {\left (b-\frac {b^2-2 a c}{\sqrt {b^2-4 a c}}\right ) \log \left (-\sqrt [3]{2} \sqrt [3]{c} x \sqrt [3]{b-\sqrt {b^2-4 a c}}+\left (b-\sqrt {b^2-4 a c}\right )^{2/3}+2^{2/3} c^{2/3} x^2\right )}{6 \sqrt [3]{2} c^{4/3} \left (b-\sqrt {b^2-4 a c}\right )^{2/3}}+\frac {\left (\frac {b^2-2 a c}{\sqrt {b^2-4 a c}}+b\right ) \log \left (-\sqrt [3]{2} \sqrt [3]{c} x \sqrt [3]{\sqrt {b^2-4 a c}+b}+\left (\sqrt {b^2-4 a c}+b\right )^{2/3}+2^{2/3} c^{2/3} x^2\right )}{6 \sqrt [3]{2} c^{4/3} \left (\sqrt {b^2-4 a c}+b\right )^{2/3}}-\frac {\left (b-\frac {b^2-2 a c}{\sqrt {b^2-4 a c}}\right ) \log \left (\sqrt [3]{b-\sqrt {b^2-4 a c}}+\sqrt [3]{2} \sqrt [3]{c} x\right )}{3 \sqrt [3]{2} c^{4/3} \left (b-\sqrt {b^2-4 a c}\right )^{2/3}}-\frac {\left (\frac {b^2-2 a c}{\sqrt {b^2-4 a c}}+b\right ) \log \left (\sqrt [3]{\sqrt {b^2-4 a c}+b}+\sqrt [3]{2} \sqrt [3]{c} x\right )}{3 \sqrt [3]{2} c^{4/3} \left (\sqrt {b^2-4 a c}+b\right )^{2/3}}+\frac {x}{c} \end {gather*}

Antiderivative was successfully verified.

[In]

Int[(c + a/x^6 + b/x^3)^(-1),x]

[Out]

x/c + ((b - (b^2 - 2*a*c)/Sqrt[b^2 - 4*a*c])*ArcTan[(1 - (2*2^(1/3)*c^(1/3)*x)/(b - Sqrt[b^2 - 4*a*c])^(1/3))/
Sqrt[3]])/(2^(1/3)*Sqrt[3]*c^(4/3)*(b - Sqrt[b^2 - 4*a*c])^(2/3)) + ((b + (b^2 - 2*a*c)/Sqrt[b^2 - 4*a*c])*Arc
Tan[(1 - (2*2^(1/3)*c^(1/3)*x)/(b + Sqrt[b^2 - 4*a*c])^(1/3))/Sqrt[3]])/(2^(1/3)*Sqrt[3]*c^(4/3)*(b + Sqrt[b^2
- 4*a*c])^(2/3)) - ((b - (b^2 - 2*a*c)/Sqrt[b^2 - 4*a*c])*Log[(b - Sqrt[b^2 - 4*a*c])^(1/3) + 2^(1/3)*c^(1/3)
*x])/(3*2^(1/3)*c^(4/3)*(b - Sqrt[b^2 - 4*a*c])^(2/3)) - ((b + (b^2 - 2*a*c)/Sqrt[b^2 - 4*a*c])*Log[(b + Sqrt[
b^2 - 4*a*c])^(1/3) + 2^(1/3)*c^(1/3)*x])/(3*2^(1/3)*c^(4/3)*(b + Sqrt[b^2 - 4*a*c])^(2/3)) + ((b - (b^2 - 2*a
*c)/Sqrt[b^2 - 4*a*c])*Log[(b - Sqrt[b^2 - 4*a*c])^(2/3) - 2^(1/3)*c^(1/3)*(b - Sqrt[b^2 - 4*a*c])^(1/3)*x + 2
^(2/3)*c^(2/3)*x^2])/(6*2^(1/3)*c^(4/3)*(b - Sqrt[b^2 - 4*a*c])^(2/3)) + ((b + (b^2 - 2*a*c)/Sqrt[b^2 - 4*a*c]
)*Log[(b + Sqrt[b^2 - 4*a*c])^(2/3) - 2^(1/3)*c^(1/3)*(b + Sqrt[b^2 - 4*a*c])^(1/3)*x + 2^(2/3)*c^(2/3)*x^2])/
(6*2^(1/3)*c^(4/3)*(b + Sqrt[b^2 - 4*a*c])^(2/3))

Rule 31

Int[((a_) + (b_.)*(x_))^(-1), x_Symbol] :> Simp[Log[RemoveContent[a + b*x, x]]/b, x] /; FreeQ[{a, b}, x]

Rule 206

Int[((a_) + (b_.)*(x_)^3)^(-1), x_Symbol] :> Dist[1/(3*Rt[a, 3]^2), Int[1/(Rt[a, 3] + Rt[b, 3]*x), x], x] + Di
st[1/(3*Rt[a, 3]^2), Int[(2*Rt[a, 3] - Rt[b, 3]*x)/(Rt[a, 3]^2 - Rt[a, 3]*Rt[b, 3]*x + Rt[b, 3]^2*x^2), x], x]
/; FreeQ[{a, b}, x]

Rule 210

Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(-(Rt[-a, 2]*Rt[-b, 2])^(-1))*ArcTan[Rt[-b, 2]*(x/Rt[-a, 2])
], x] /; FreeQ[{a, b}, x] && PosQ[a/b] && (LtQ[a, 0] || LtQ[b, 0])

Rule 631

Int[((a_) + (b_.)*(x_) + (c_.)*(x_)^2)^(-1), x_Symbol] :> With[{q = 1 - 4*Simplify[a*(c/b^2)]}, Dist[-2/b, Sub
st[Int[1/(q - x^2), x], x, 1 + 2*c*(x/b)], x] /; RationalQ[q] && (EqQ[q^2, 1] ||  !RationalQ[b^2 - 4*a*c])] /;
FreeQ[{a, b, c}, x] && NeQ[b^2 - 4*a*c, 0]

Rule 642

Int[((d_) + (e_.)*(x_))/((a_.) + (b_.)*(x_) + (c_.)*(x_)^2), x_Symbol] :> Simp[d*(Log[RemoveContent[a + b*x +
c*x^2, x]]/b), x] /; FreeQ[{a, b, c, d, e}, x] && EqQ[2*c*d - b*e, 0]

Rule 648

Int[((d_.) + (e_.)*(x_))/((a_) + (b_.)*(x_) + (c_.)*(x_)^2), x_Symbol] :> Dist[(2*c*d - b*e)/(2*c), Int[1/(a +
b*x + c*x^2), x], x] + Dist[e/(2*c), Int[(b + 2*c*x)/(a + b*x + c*x^2), x], x] /; FreeQ[{a, b, c, d, e}, x] &
& NeQ[2*c*d - b*e, 0] && NeQ[b^2 - 4*a*c, 0] &&  !NiceSqrtQ[b^2 - 4*a*c]

Rule 1354

Int[((a_) + (c_.)*(x_)^(n2_.) + (b_.)*(x_)^(n_))^(p_.), x_Symbol] :> Int[x^(2*n*p)*(c + b/x^n + a/x^(2*n))^p,
x] /; FreeQ[{a, b, c}, x] && EqQ[n2, 2*n] && LtQ[n, 0] && IntegerQ[p]

Rule 1381

Int[((d_.)*(x_))^(m_.)*((a_) + (c_.)*(x_)^(n2_.) + (b_.)*(x_)^(n_))^(p_), x_Symbol] :> Simp[d^(2*n - 1)*(d*x)^
(m - 2*n + 1)*((a + b*x^n + c*x^(2*n))^(p + 1)/(c*(m + 2*n*p + 1))), x] - Dist[d^(2*n)/(c*(m + 2*n*p + 1)), In
t[(d*x)^(m - 2*n)*Simp[a*(m - 2*n + 1) + b*(m + n*(p - 1) + 1)*x^n, x]*(a + b*x^n + c*x^(2*n))^p, x], x] /; Fr
eeQ[{a, b, c, d, p}, x] && EqQ[n2, 2*n] && NeQ[b^2 - 4*a*c, 0] && IGtQ[n, 0] && GtQ[m, 2*n - 1] && NeQ[m + 2*n
*p + 1, 0] && IntegerQ[p]

Rule 1436

Int[((d_) + (e_.)*(x_)^(n_))/((a_) + (b_.)*(x_)^(n_) + (c_.)*(x_)^(n2_)), x_Symbol] :> With[{q = Rt[b^2 - 4*a*
c, 2]}, Dist[e/2 + (2*c*d - b*e)/(2*q), Int[1/(b/2 - q/2 + c*x^n), x], x] + Dist[e/2 - (2*c*d - b*e)/(2*q), In
t[1/(b/2 + q/2 + c*x^n), x], x]] /; FreeQ[{a, b, c, d, e, n}, x] && EqQ[n2, 2*n] && NeQ[b^2 - 4*a*c, 0] && NeQ
[c*d^2 - b*d*e + a*e^2, 0] && (PosQ[b^2 - 4*a*c] ||  !IGtQ[n/2, 0])

Rubi steps

\begin {align*} \int \frac {1}{c+\frac {a}{x^6}+\frac {b}{x^3}} \, dx &=\int \frac {x^6}{a+b x^3+c x^6} \, dx\\ &=\frac {x}{c}-\frac {\int \frac {a+b x^3}{a+b x^3+c x^6} \, dx}{c}\\ &=\frac {x}{c}-\frac {\left (b-\frac {b^2-2 a c}{\sqrt {b^2-4 a c}}\right ) \int \frac {1}{\frac {b}{2}-\frac {1}{2} \sqrt {b^2-4 a c}+c x^3} \, dx}{2 c}-\frac {\left (b+\frac {b^2-2 a c}{\sqrt {b^2-4 a c}}\right ) \int \frac {1}{\frac {b}{2}+\frac {1}{2} \sqrt {b^2-4 a c}+c x^3} \, dx}{2 c}\\ &=\frac {x}{c}-\frac {\left (b-\frac {b^2-2 a c}{\sqrt {b^2-4 a c}}\right ) \int \frac {1}{\frac {\sqrt [3]{b-\sqrt {b^2-4 a c}}}{\sqrt [3]{2}}+\sqrt [3]{c} x} \, dx}{3 \sqrt [3]{2} c \left (b-\sqrt {b^2-4 a c}\right )^{2/3}}-\frac {\left (b-\frac {b^2-2 a c}{\sqrt {b^2-4 a c}}\right ) \int \frac {2^{2/3} \sqrt [3]{b-\sqrt {b^2-4 a c}}-\sqrt [3]{c} x}{\frac {\left (b-\sqrt {b^2-4 a c}\right )^{2/3}}{2^{2/3}}-\frac {\sqrt [3]{c} \sqrt [3]{b-\sqrt {b^2-4 a c}} x}{\sqrt [3]{2}}+c^{2/3} x^2} \, dx}{3 \sqrt [3]{2} c \left (b-\sqrt {b^2-4 a c}\right )^{2/3}}-\frac {\left (b+\frac {b^2-2 a c}{\sqrt {b^2-4 a c}}\right ) \int \frac {1}{\frac {\sqrt [3]{b+\sqrt {b^2-4 a c}}}{\sqrt [3]{2}}+\sqrt [3]{c} x} \, dx}{3 \sqrt [3]{2} c \left (b+\sqrt {b^2-4 a c}\right )^{2/3}}-\frac {\left (b+\frac {b^2-2 a c}{\sqrt {b^2-4 a c}}\right ) \int \frac {2^{2/3} \sqrt [3]{b+\sqrt {b^2-4 a c}}-\sqrt [3]{c} x}{\frac {\left (b+\sqrt {b^2-4 a c}\right )^{2/3}}{2^{2/3}}-\frac {\sqrt [3]{c} \sqrt [3]{b+\sqrt {b^2-4 a c}} x}{\sqrt [3]{2}}+c^{2/3} x^2} \, dx}{3 \sqrt [3]{2} c \left (b+\sqrt {b^2-4 a c}\right )^{2/3}}\\ &=\frac {x}{c}-\frac {\left (b-\frac {b^2-2 a c}{\sqrt {b^2-4 a c}}\right ) \log \left (\sqrt [3]{b-\sqrt {b^2-4 a c}}+\sqrt [3]{2} \sqrt [3]{c} x\right )}{3 \sqrt [3]{2} c^{4/3} \left (b-\sqrt {b^2-4 a c}\right )^{2/3}}-\frac {\left (b+\frac {b^2-2 a c}{\sqrt {b^2-4 a c}}\right ) \log \left (\sqrt [3]{b+\sqrt {b^2-4 a c}}+\sqrt [3]{2} \sqrt [3]{c} x\right )}{3 \sqrt [3]{2} c^{4/3} \left (b+\sqrt {b^2-4 a c}\right )^{2/3}}+\frac {\left (b-\frac {b^2-2 a c}{\sqrt {b^2-4 a c}}\right ) \int \frac {-\frac {\sqrt [3]{c} \sqrt [3]{b-\sqrt {b^2-4 a c}}}{\sqrt [3]{2}}+2 c^{2/3} x}{\frac {\left (b-\sqrt {b^2-4 a c}\right )^{2/3}}{2^{2/3}}-\frac {\sqrt [3]{c} \sqrt [3]{b-\sqrt {b^2-4 a c}} x}{\sqrt [3]{2}}+c^{2/3} x^2} \, dx}{6 \sqrt [3]{2} c^{4/3} \left (b-\sqrt {b^2-4 a c}\right )^{2/3}}-\frac {\left (b-\frac {b^2-2 a c}{\sqrt {b^2-4 a c}}\right ) \int \frac {1}{\frac {\left (b-\sqrt {b^2-4 a c}\right )^{2/3}}{2^{2/3}}-\frac {\sqrt [3]{c} \sqrt [3]{b-\sqrt {b^2-4 a c}} x}{\sqrt [3]{2}}+c^{2/3} x^2} \, dx}{2\ 2^{2/3} c \sqrt [3]{b-\sqrt {b^2-4 a c}}}+\frac {\left (b+\frac {b^2-2 a c}{\sqrt {b^2-4 a c}}\right ) \int \frac {-\frac {\sqrt [3]{c} \sqrt [3]{b+\sqrt {b^2-4 a c}}}{\sqrt [3]{2}}+2 c^{2/3} x}{\frac {\left (b+\sqrt {b^2-4 a c}\right )^{2/3}}{2^{2/3}}-\frac {\sqrt [3]{c} \sqrt [3]{b+\sqrt {b^2-4 a c}} x}{\sqrt [3]{2}}+c^{2/3} x^2} \, dx}{6 \sqrt [3]{2} c^{4/3} \left (b+\sqrt {b^2-4 a c}\right )^{2/3}}-\frac {\left (b+\frac {b^2-2 a c}{\sqrt {b^2-4 a c}}\right ) \int \frac {1}{\frac {\left (b+\sqrt {b^2-4 a c}\right )^{2/3}}{2^{2/3}}-\frac {\sqrt [3]{c} \sqrt [3]{b+\sqrt {b^2-4 a c}} x}{\sqrt [3]{2}}+c^{2/3} x^2} \, dx}{2\ 2^{2/3} c \sqrt [3]{b+\sqrt {b^2-4 a c}}}\\ &=\frac {x}{c}-\frac {\left (b-\frac {b^2-2 a c}{\sqrt {b^2-4 a c}}\right ) \log \left (\sqrt [3]{b-\sqrt {b^2-4 a c}}+\sqrt [3]{2} \sqrt [3]{c} x\right )}{3 \sqrt [3]{2} c^{4/3} \left (b-\sqrt {b^2-4 a c}\right )^{2/3}}-\frac {\left (b+\frac {b^2-2 a c}{\sqrt {b^2-4 a c}}\right ) \log \left (\sqrt [3]{b+\sqrt {b^2-4 a c}}+\sqrt [3]{2} \sqrt [3]{c} x\right )}{3 \sqrt [3]{2} c^{4/3} \left (b+\sqrt {b^2-4 a c}\right )^{2/3}}+\frac {\left (b-\frac {b^2-2 a c}{\sqrt {b^2-4 a c}}\right ) \log \left (\left (b-\sqrt {b^2-4 a c}\right )^{2/3}-\sqrt [3]{2} \sqrt [3]{c} \sqrt [3]{b-\sqrt {b^2-4 a c}} x+2^{2/3} c^{2/3} x^2\right )}{6 \sqrt [3]{2} c^{4/3} \left (b-\sqrt {b^2-4 a c}\right )^{2/3}}+\frac {\left (b+\frac {b^2-2 a c}{\sqrt {b^2-4 a c}}\right ) \log \left (\left (b+\sqrt {b^2-4 a c}\right )^{2/3}-\sqrt [3]{2} \sqrt [3]{c} \sqrt [3]{b+\sqrt {b^2-4 a c}} x+2^{2/3} c^{2/3} x^2\right )}{6 \sqrt [3]{2} c^{4/3} \left (b+\sqrt {b^2-4 a c}\right )^{2/3}}-\frac {\left (b-\frac {b^2-2 a c}{\sqrt {b^2-4 a c}}\right ) \text {Subst}\left (\int \frac {1}{-3-x^2} \, dx,x,1-\frac {2 \sqrt [3]{2} \sqrt [3]{c} x}{\sqrt [3]{b-\sqrt {b^2-4 a c}}}\right )}{\sqrt [3]{2} c^{4/3} \left (b-\sqrt {b^2-4 a c}\right )^{2/3}}-\frac {\left (b+\frac {b^2-2 a c}{\sqrt {b^2-4 a c}}\right ) \text {Subst}\left (\int \frac {1}{-3-x^2} \, dx,x,1-\frac {2 \sqrt [3]{2} \sqrt [3]{c} x}{\sqrt [3]{b+\sqrt {b^2-4 a c}}}\right )}{\sqrt [3]{2} c^{4/3} \left (b+\sqrt {b^2-4 a c}\right )^{2/3}}\\ &=\frac {x}{c}+\frac {\left (b-\frac {b^2-2 a c}{\sqrt {b^2-4 a c}}\right ) \tan ^{-1}\left (\frac {1-\frac {2 \sqrt [3]{2} \sqrt [3]{c} x}{\sqrt [3]{b-\sqrt {b^2-4 a c}}}}{\sqrt {3}}\right )}{\sqrt [3]{2} \sqrt {3} c^{4/3} \left (b-\sqrt {b^2-4 a c}\right )^{2/3}}+\frac {\left (b+\frac {b^2-2 a c}{\sqrt {b^2-4 a c}}\right ) \tan ^{-1}\left (\frac {1-\frac {2 \sqrt [3]{2} \sqrt [3]{c} x}{\sqrt [3]{b+\sqrt {b^2-4 a c}}}}{\sqrt {3}}\right )}{\sqrt [3]{2} \sqrt {3} c^{4/3} \left (b+\sqrt {b^2-4 a c}\right )^{2/3}}-\frac {\left (b-\frac {b^2-2 a c}{\sqrt {b^2-4 a c}}\right ) \log \left (\sqrt [3]{b-\sqrt {b^2-4 a c}}+\sqrt [3]{2} \sqrt [3]{c} x\right )}{3 \sqrt [3]{2} c^{4/3} \left (b-\sqrt {b^2-4 a c}\right )^{2/3}}-\frac {\left (b+\frac {b^2-2 a c}{\sqrt {b^2-4 a c}}\right ) \log \left (\sqrt [3]{b+\sqrt {b^2-4 a c}}+\sqrt [3]{2} \sqrt [3]{c} x\right )}{3 \sqrt [3]{2} c^{4/3} \left (b+\sqrt {b^2-4 a c}\right )^{2/3}}+\frac {\left (b-\frac {b^2-2 a c}{\sqrt {b^2-4 a c}}\right ) \log \left (\left (b-\sqrt {b^2-4 a c}\right )^{2/3}-\sqrt [3]{2} \sqrt [3]{c} \sqrt [3]{b-\sqrt {b^2-4 a c}} x+2^{2/3} c^{2/3} x^2\right )}{6 \sqrt [3]{2} c^{4/3} \left (b-\sqrt {b^2-4 a c}\right )^{2/3}}+\frac {\left (b+\frac {b^2-2 a c}{\sqrt {b^2-4 a c}}\right ) \log \left (\left (b+\sqrt {b^2-4 a c}\right )^{2/3}-\sqrt [3]{2} \sqrt [3]{c} \sqrt [3]{b+\sqrt {b^2-4 a c}} x+2^{2/3} c^{2/3} x^2\right )}{6 \sqrt [3]{2} c^{4/3} \left (b+\sqrt {b^2-4 a c}\right )^{2/3}}\\ \end {align*}

________________________________________________________________________________________

Mathematica [C] Result contains higher order function than in optimal. Order 9 vs. order 3 in optimal.
time = 0.02, size = 70, normalized size = 0.11 \begin {gather*} \frac {x}{c}-\frac {\text {RootSum}\left [a+b \text {$\#$1}^3+c \text {$\#$1}^6\&,\frac {a \log (x-\text {$\#$1})+b \log (x-\text {$\#$1}) \text {$\#$1}^3}{b \text {$\#$1}^2+2 c \text {$\#$1}^5}\&\right ]}{3 c} \end {gather*}

Antiderivative was successfully verified.

[In]

Integrate[(c + a/x^6 + b/x^3)^(-1),x]

[Out]

x/c - RootSum[a + b*#1^3 + c*#1^6 & , (a*Log[x - #1] + b*Log[x - #1]*#1^3)/(b*#1^2 + 2*c*#1^5) & ]/(3*c)

________________________________________________________________________________________

Maple [C] Result contains higher order function than in optimal. Order 9 vs. order 3.
time = 0.20, size = 59, normalized size = 0.09

method	result	size
default	\(\frac {x}{c}+\frac {\munderset {\textit {\_R} =\RootOf \left (\textit {\_Z}^{6} c +\textit {\_Z}^{3} b +a \right )}{\sum }\frac {\left (-\textit {\_R}^{3} b -a \right ) \ln \left (x -\textit {\_R} \right )}{2 \textit {\_R}^{5} c +b \,\textit {\_R}^{2}}}{3 c}\)	\(59\)
risch	\(\frac {x}{c}+\frac {\munderset {\textit {\_R} =\RootOf \left (\textit {\_Z}^{6} c +\textit {\_Z}^{3} b +a \right )}{\sum }\frac {\left (-\textit {\_R}^{3} b -a \right ) \ln \left (x -\textit {\_R} \right )}{2 \textit {\_R}^{5} c +b \,\textit {\_R}^{2}}}{3 c}\)	\(59\)

Verification of antiderivative is not currently implemented for this CAS.

[In]

int(1/(c+a/x^6+b/x^3),x,method=_RETURNVERBOSE)

[Out]

x/c+1/3/c*sum((-_R^3*b-a)/(2*_R^5*c+_R^2*b)*ln(x-_R),_R=RootOf(_Z^6*c+_Z^3*b+a))

________________________________________________________________________________________

Maxima [F]
time = 0.00, size = 0, normalized size = 0.00 \begin {gather*} \text {Failed to integrate} \end {gather*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(1/(c+a/x^6+b/x^3),x, algorithm="maxima")

[Out]

x/c - integrate((b*x^3 + a)/(c*x^6 + b*x^3 + a), x)/c

________________________________________________________________________________________

Fricas [B] Leaf count of result is larger than twice the leaf count of optimal. 5841 vs. \(2 (495) = 990\).
time = 1.67, size = 5841, normalized size = 9.26 \begin {gather*} \text {Too large to display} \end {gather*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(1/(c+a/x^6+b/x^3),x, algorithm="fricas")

[Out]

1/6*(4*sqrt(3)*(1/2)^(1/3)*c*(-(b^3 - 2*a*b*c + (b^2*c^4 - 4*a*c^5)*sqrt((b^8 - 8*a*b^6*c + 20*a^2*b^4*c^2 - 1
6*a^3*b^2*c^3 + 4*a^4*c^4)/(b^6*c^8 - 12*a*b^4*c^9 + 48*a^2*b^2*c^10 - 64*a^3*c^11)))/(b^2*c^4 - 4*a*c^5))^(1/
3)*arctan(1/6*(sqrt(2)*(1/2)^(2/3)*sqrt(2*(a^2*b^8 - 8*a^3*b^6*c + 20*a^4*b^4*c^2 - 16*a^5*b^2*c^3 + 4*a^6*c^4
)*x^2 + (1/2)^(2/3)*(b^12 - 14*a*b^10*c + 76*a^2*b^8*c^2 - 200*a^3*b^6*c^3 + 260*a^4*b^4*c^4 - 152*a^5*b^2*c^5
+ 32*a^6*c^6 - (b^11*c^4 - 16*a*b^9*c^5 + 98*a^2*b^7*c^6 - 280*a^3*b^5*c^7 + 352*a^4*b^3*c^8 - 128*a^5*b*c^9)
*sqrt((b^8 - 8*a*b^6*c + 20*a^2*b^4*c^2 - 16*a^3*b^2*c^3 + 4*a^4*c^4)/(b^6*c^8 - 12*a*b^4*c^9 + 48*a^2*b^2*c^1
0 - 64*a^3*c^11)))*(-(b^3 - 2*a*b*c + (b^2*c^4 - 4*a*c^5)*sqrt((b^8 - 8*a*b^6*c + 20*a^2*b^4*c^2 - 16*a^3*b^2*
c^3 + 4*a^4*c^4)/(b^6*c^8 - 12*a*b^4*c^9 + 48*a^2*b^2*c^10 - 64*a^3*c^11)))/(b^2*c^4 - 4*a*c^5))^(2/3) + (1/2)
^(1/3)*((a*b^9*c^4 - 12*a^2*b^7*c^5 + 50*a^3*b^5*c^6 - 80*a^4*b^3*c^7 + 32*a^5*b*c^8)*x*sqrt((b^8 - 8*a*b^6*c
+ 20*a^2*b^4*c^2 - 16*a^3*b^2*c^3 + 4*a^4*c^4)/(b^6*c^8 - 12*a*b^4*c^9 + 48*a^2*b^2*c^10 - 64*a^3*c^11)) - (a*
b^10 - 12*a^2*b^8*c + 52*a^3*b^6*c^2 - 96*a^4*b^4*c^3 + 68*a^5*b^2*c^4 - 16*a^6*c^5)*x)*(-(b^3 - 2*a*b*c + (b^
2*c^4 - 4*a*c^5)*sqrt((b^8 - 8*a*b^6*c + 20*a^2*b^4*c^2 - 16*a^3*b^2*c^3 + 4*a^4*c^4)/(b^6*c^8 - 12*a*b^4*c^9
+ 48*a^2*b^2*c^10 - 64*a^3*c^11)))/(b^2*c^4 - 4*a*c^5))^(1/3))*(sqrt(3)*(b^8*c^4 - 13*a*b^6*c^5 + 60*a^2*b^4*c
^6 - 112*a^3*b^2*c^7 + 64*a^4*c^8)*sqrt((b^8 - 8*a*b^6*c + 20*a^2*b^4*c^2 - 16*a^3*b^2*c^3 + 4*a^4*c^4)/(b^6*c
^8 - 12*a*b^4*c^9 + 48*a^2*b^2*c^10 - 64*a^3*c^11)) - sqrt(3)*(b^9 - 11*a*b^7*c + 42*a^2*b^5*c^2 - 62*a^3*b^3*
c^3 + 24*a^4*b*c^4))*(-(b^3 - 2*a*b*c + (b^2*c^4 - 4*a*c^5)*sqrt((b^8 - 8*a*b^6*c + 20*a^2*b^4*c^2 - 16*a^3*b^
2*c^3 + 4*a^4*c^4)/(b^6*c^8 - 12*a*b^4*c^9 + 48*a^2*b^2*c^10 - 64*a^3*c^11)))/(b^2*c^4 - 4*a*c^5))^(2/3) - 2*(
1/2)^(2/3)*(sqrt(3)*(a*b^12*c^4 - 17*a^2*b^10*c^5 + 114*a^3*b^8*c^6 - 378*a^4*b^6*c^7 + 632*a^5*b^4*c^8 - 480*
a^6*b^2*c^9 + 128*a^7*c^10)*x*sqrt((b^8 - 8*a*b^6*c + 20*a^2*b^4*c^2 - 16*a^3*b^2*c^3 + 4*a^4*c^4)/(b^6*c^8 -
12*a*b^4*c^9 + 48*a^2*b^2*c^10 - 64*a^3*c^11)) - sqrt(3)*(a*b^13 - 15*a^2*b^11*c + 88*a^3*b^9*c^2 - 252*a^4*b^
7*c^3 + 356*a^5*b^5*c^4 - 220*a^6*b^3*c^5 + 48*a^7*b*c^6)*x)*(-(b^3 - 2*a*b*c + (b^2*c^4 - 4*a*c^5)*sqrt((b^8
- 8*a*b^6*c + 20*a^2*b^4*c^2 - 16*a^3*b^2*c^3 + 4*a^4*c^4)/(b^6*c^8 - 12*a*b^4*c^9 + 48*a^2*b^2*c^10 - 64*a^3*
c^11)))/(b^2*c^4 - 4*a*c^5))^(2/3) - 2*sqrt(3)*(a^4*b^8 - 8*a^5*b^6*c + 20*a^6*b^4*c^2 - 16*a^7*b^2*c^3 + 4*a^
8*c^4))/(a^4*b^8 - 8*a^5*b^6*c + 20*a^6*b^4*c^2 - 16*a^7*b^2*c^3 + 4*a^8*c^4)) - 4*sqrt(3)*(1/2)^(1/3)*c*(-(b^
3 - 2*a*b*c - (b^2*c^4 - 4*a*c^5)*sqrt((b^8 - 8*a*b^6*c + 20*a^2*b^4*c^2 - 16*a^3*b^2*c^3 + 4*a^4*c^4)/(b^6*c^
8 - 12*a*b^4*c^9 + 48*a^2*b^2*c^10 - 64*a^3*c^11)))/(b^2*c^4 - 4*a*c^5))^(1/3)*arctan(1/6*(sqrt(2)*(1/2)^(2/3)
*sqrt(2*(a^2*b^8 - 8*a^3*b^6*c + 20*a^4*b^4*c^2 - 16*a^5*b^2*c^3 + 4*a^6*c^4)*x^2 + (1/2)^(2/3)*(b^12 - 14*a*b
^10*c + 76*a^2*b^8*c^2 - 200*a^3*b^6*c^3 + 260*a^4*b^4*c^4 - 152*a^5*b^2*c^5 + 32*a^6*c^6 + (b^11*c^4 - 16*a*b
^9*c^5 + 98*a^2*b^7*c^6 - 280*a^3*b^5*c^7 + 352*a^4*b^3*c^8 - 128*a^5*b*c^9)*sqrt((b^8 - 8*a*b^6*c + 20*a^2*b^
4*c^2 - 16*a^3*b^2*c^3 + 4*a^4*c^4)/(b^6*c^8 - 12*a*b^4*c^9 + 48*a^2*b^2*c^10 - 64*a^3*c^11)))*(-(b^3 - 2*a*b*
c - (b^2*c^4 - 4*a*c^5)*sqrt((b^8 - 8*a*b^6*c + 20*a^2*b^4*c^2 - 16*a^3*b^2*c^3 + 4*a^4*c^4)/(b^6*c^8 - 12*a*b
^4*c^9 + 48*a^2*b^2*c^10 - 64*a^3*c^11)))/(b^2*c^4 - 4*a*c^5))^(2/3) - (1/2)^(1/3)*((a*b^9*c^4 - 12*a^2*b^7*c^
5 + 50*a^3*b^5*c^6 - 80*a^4*b^3*c^7 + 32*a^5*b*c^8)*x*sqrt((b^8 - 8*a*b^6*c + 20*a^2*b^4*c^2 - 16*a^3*b^2*c^3
+ 4*a^4*c^4)/(b^6*c^8 - 12*a*b^4*c^9 + 48*a^2*b^2*c^10 - 64*a^3*c^11)) + (a*b^10 - 12*a^2*b^8*c + 52*a^3*b^6*c
^2 - 96*a^4*b^4*c^3 + 68*a^5*b^2*c^4 - 16*a^6*c^5)*x)*(-(b^3 - 2*a*b*c - (b^2*c^4 - 4*a*c^5)*sqrt((b^8 - 8*a*b
^6*c + 20*a^2*b^4*c^2 - 16*a^3*b^2*c^3 + 4*a^4*c^4)/(b^6*c^8 - 12*a*b^4*c^9 + 48*a^2*b^2*c^10 - 64*a^3*c^11)))
/(b^2*c^4 - 4*a*c^5))^(1/3))*(sqrt(3)*(b^8*c^4 - 13*a*b^6*c^5 + 60*a^2*b^4*c^6 - 112*a^3*b^2*c^7 + 64*a^4*c^8)
*sqrt((b^8 - 8*a*b^6*c + 20*a^2*b^4*c^2 - 16*a^3*b^2*c^3 + 4*a^4*c^4)/(b^6*c^8 - 12*a*b^4*c^9 + 48*a^2*b^2*c^1
0 - 64*a^3*c^11)) + sqrt(3)*(b^9 - 11*a*b^7*c + 42*a^2*b^5*c^2 - 62*a^3*b^3*c^3 + 24*a^4*b*c^4))*(-(b^3 - 2*a*
b*c - (b^2*c^4 - 4*a*c^5)*sqrt((b^8 - 8*a*b^6*c + 20*a^2*b^4*c^2 - 16*a^3*b^2*c^3 + 4*a^4*c^4)/(b^6*c^8 - 12*a
*b^4*c^9 + 48*a^2*b^2*c^10 - 64*a^3*c^11)))/(b^2*c^4 - 4*a*c^5))^(2/3) - 2*(1/2)^(2/3)*(sqrt(3)*(a*b^12*c^4 -
17*a^2*b^10*c^5 + 114*a^3*b^8*c^6 - 378*a^4*b^6*c^7 + 632*a^5*b^4*c^8 - 480*a^6*b^2*c^9 + 128*a^7*c^10)*x*sqrt
((b^8 - 8*a*b^6*c + 20*a^2*b^4*c^2 - 16*a^3*b^2*c^3 + 4*a^4*c^4)/(b^6*c^8 - 12*a*b^4*c^9 + 48*a^2*b^2*c^10 - 6
4*a^3*c^11)) + sqrt(3)*(a*b^13 - 15*a^2*b^11*c + 88*a^3*b^9*c^2 - 252*a^4*b^7*c^3 + 356*a^5*b^5*c^4 - 220*a^6*
b^3*c^5 + 48*a^7*b*c^6)*x)*(-(b^3 - 2*a*b*c - (b^2*c^4 - 4*a*c^5)*sqrt((b^8 - 8*a*b^6*c + 20*a^2*b^4*c^2 - 16*
a^3*b^2*c^3 + 4*a^4*c^4)/(b^6*c^8 - 12*a*b^4*c^...

________________________________________________________________________________________

Sympy [A]
time = 52.83, size = 196, normalized size = 0.31 \begin {gather*} \operatorname {RootSum} {\left (t^{6} \cdot \left (46656 a^{3} c^{7} - 34992 a^{2} b^{2} c^{6} + 8748 a b^{4} c^{5} - 729 b^{6} c^{4}\right ) + t^{3} \cdot \left (864 a^{3} b c^{3} - 864 a^{2} b^{3} c^{2} + 270 a b^{5} c - 27 b^{7}\right ) + a^{4}, \left ( t \mapsto t \log {\left (x + \frac {1296 t^{4} a^{2} b c^{6} - 648 t^{4} a b^{3} c^{5} + 81 t^{4} b^{5} c^{4} - 12 t a^{3} c^{3} + 39 t a^{2} b^{2} c^{2} - 21 t a b^{4} c + 3 t b^{6}}{2 a^{3} c^{2} - 4 a^{2} b^{2} c + a b^{4}} \right )} \right )\right )} + \frac {x}{c} \end {gather*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(1/(c+a/x**6+b/x**3),x)

[Out]

RootSum(_t**6*(46656*a**3*c**7 - 34992*a**2*b**2*c**6 + 8748*a*b**4*c**5 - 729*b**6*c**4) + _t**3*(864*a**3*b*
c**3 - 864*a**2*b**3*c**2 + 270*a*b**5*c - 27*b**7) + a**4, Lambda(_t, _t*log(x + (1296*_t**4*a**2*b*c**6 - 64
8*_t**4*a*b**3*c**5 + 81*_t**4*b**5*c**4 - 12*_t*a**3*c**3 + 39*_t*a**2*b**2*c**2 - 21*_t*a*b**4*c + 3*_t*b**6
)/(2*a**3*c**2 - 4*a**2*b**2*c + a*b**4)))) + x/c

________________________________________________________________________________________

Giac [F]
time = 0.00, size = 0, normalized size = 0.00 \begin {gather*} \text {could not integrate} \end {gather*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(1/(c+a/x^6+b/x^3),x, algorithm="giac")

[Out]

integrate(1/(c + b/x^3 + a/x^6), x)

________________________________________________________________________________________

Mupad [B]
time = 4.54, size = 2280, normalized size = 3.61 \begin {gather*} \text {Too large to display} \end {gather*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

int(1/(c + a/x^6 + b/x^3),x)

[Out]

log((3*a^2*x*(b^4 + 2*a^2*c^2 - 4*a*b^2*c))/c - (3*2^(2/3)*a*(-(b^4*(-(4*a*c - b^2)^3)^(1/2) - b^7 + 32*a^3*b*
c^3 - 32*a^2*b^3*c^2 + 2*a^2*c^2*(-(4*a*c - b^2)^3)^(1/2) + 10*a*b^5*c - 4*a*b^2*c*(-(4*a*c - b^2)^3)^(1/2))/(
c^4*(4*a*c - b^2)^3))^(1/3)*(b^4 + 2*a^2*c^2 - 4*a*b^2*c)*(b*(-(4*a*c - b^2)^3)^(1/2) + b^4 + 16*a^2*c^2 - 8*a
*b^2*c))/(4*c*(4*a*c - b^2)))*(-(b^4*(-(4*a*c - b^2)^3)^(1/2) - b^7 + 32*a^3*b*c^3 - 32*a^2*b^3*c^2 + 2*a^2*c^
2*(-(4*a*c - b^2)^3)^(1/2) + 10*a*b^5*c - 4*a*b^2*c*(-(4*a*c - b^2)^3)^(1/2))/(54*(64*a^3*c^7 - b^6*c^4 + 12*a
*b^4*c^5 - 48*a^2*b^2*c^6)))^(1/3) + x/c + log((3*a^2*x*(b^4 + 2*a^2*c^2 - 4*a*b^2*c))/c + (3*2^(2/3)*a*((b^7
+ b^4*(-(4*a*c - b^2)^3)^(1/2) - 32*a^3*b*c^3 + 32*a^2*b^3*c^2 + 2*a^2*c^2*(-(4*a*c - b^2)^3)^(1/2) - 10*a*b^5
*c - 4*a*b^2*c*(-(4*a*c - b^2)^3)^(1/2))/(c^4*(4*a*c - b^2)^3))^(1/3)*(b^4 + 2*a^2*c^2 - 4*a*b^2*c)*(b*(-(4*a*
c - b^2)^3)^(1/2) - b^4 - 16*a^2*c^2 + 8*a*b^2*c))/(4*c*(4*a*c - b^2)))*((b^7 + b^4*(-(4*a*c - b^2)^3)^(1/2) -
32*a^3*b*c^3 + 32*a^2*b^3*c^2 + 2*a^2*c^2*(-(4*a*c - b^2)^3)^(1/2) - 10*a*b^5*c - 4*a*b^2*c*(-(4*a*c - b^2)^3
)^(1/2))/(54*(64*a^3*c^7 - b^6*c^4 + 12*a*b^4*c^5 - 48*a^2*b^2*c^6)))^(1/3) + log((3*a^2*x*(b^4 + 2*a^2*c^2 -
4*a*b^2*c))/c + (3*2^(2/3)*a*(3^(1/2)*1i - 1)*((b^7 + b^4*(-(4*a*c - b^2)^3)^(1/2) - 32*a^3*b*c^3 + 32*a^2*b^3
*c^2 + 2*a^2*c^2*(-(4*a*c - b^2)^3)^(1/2) - 10*a*b^5*c - 4*a*b^2*c*(-(4*a*c - b^2)^3)^(1/2))/(c^4*(4*a*c - b^2
)^3))^(1/3)*(b^4 + 2*a^2*c^2 - 4*a*b^2*c)*(b*(-(4*a*c - b^2)^3)^(1/2) - b^4 - 16*a^2*c^2 + 8*a*b^2*c))/(8*c*(4
*a*c - b^2)))*((3^(1/2)*1i)/2 - 1/2)*((b^7 + b^4*(-(4*a*c - b^2)^3)^(1/2) - 32*a^3*b*c^3 + 32*a^2*b^3*c^2 + 2*
a^2*c^2*(-(4*a*c - b^2)^3)^(1/2) - 10*a*b^5*c - 4*a*b^2*c*(-(4*a*c - b^2)^3)^(1/2))/(54*(64*a^3*c^7 - b^6*c^4
+ 12*a*b^4*c^5 - 48*a^2*b^2*c^6)))^(1/3) - log((3*a^2*x*(b^4 + 2*a^2*c^2 - 4*a*b^2*c))/c - (3*2^(2/3)*a*(3^(1/
2)*1i + 1)*((b^7 + b^4*(-(4*a*c - b^2)^3)^(1/2) - 32*a^3*b*c^3 + 32*a^2*b^3*c^2 + 2*a^2*c^2*(-(4*a*c - b^2)^3)
^(1/2) - 10*a*b^5*c - 4*a*b^2*c*(-(4*a*c - b^2)^3)^(1/2))/(c^4*(4*a*c - b^2)^3))^(1/3)*(b^4 + 2*a^2*c^2 - 4*a*
b^2*c)*(b*(-(4*a*c - b^2)^3)^(1/2) - b^4 - 16*a^2*c^2 + 8*a*b^2*c))/(8*c*(4*a*c - b^2)))*((3^(1/2)*1i)/2 + 1/2
)*((b^7 + b^4*(-(4*a*c - b^2)^3)^(1/2) - 32*a^3*b*c^3 + 32*a^2*b^3*c^2 + 2*a^2*c^2*(-(4*a*c - b^2)^3)^(1/2) -
10*a*b^5*c - 4*a*b^2*c*(-(4*a*c - b^2)^3)^(1/2))/(54*(64*a^3*c^7 - b^6*c^4 + 12*a*b^4*c^5 - 48*a^2*b^2*c^6)))^
(1/3) + log((3*a^2*x*(b^4 + 2*a^2*c^2 - 4*a*b^2*c))/c - (3*2^(2/3)*a*(3^(1/2)*1i - 1)*(-(b^4*(-(4*a*c - b^2)^3
)^(1/2) - b^7 + 32*a^3*b*c^3 - 32*a^2*b^3*c^2 + 2*a^2*c^2*(-(4*a*c - b^2)^3)^(1/2) + 10*a*b^5*c - 4*a*b^2*c*(-
(4*a*c - b^2)^3)^(1/2))/(c^4*(4*a*c - b^2)^3))^(1/3)*(b^4 + 2*a^2*c^2 - 4*a*b^2*c)*(b*(-(4*a*c - b^2)^3)^(1/2)
+ b^4 + 16*a^2*c^2 - 8*a*b^2*c))/(8*c*(4*a*c - b^2)))*((3^(1/2)*1i)/2 - 1/2)*(-(b^4*(-(4*a*c - b^2)^3)^(1/2)
- b^7 + 32*a^3*b*c^3 - 32*a^2*b^3*c^2 + 2*a^2*c^2*(-(4*a*c - b^2)^3)^(1/2) + 10*a*b^5*c - 4*a*b^2*c*(-(4*a*c -
b^2)^3)^(1/2))/(54*(64*a^3*c^7 - b^6*c^4 + 12*a*b^4*c^5 - 48*a^2*b^2*c^6)))^(1/3) - log((3*a^2*x*(b^4 + 2*a^2
*c^2 - 4*a*b^2*c))/c + (3*2^(2/3)*a*(3^(1/2)*1i + 1)*(-(b^4*(-(4*a*c - b^2)^3)^(1/2) - b^7 + 32*a^3*b*c^3 - 32
*a^2*b^3*c^2 + 2*a^2*c^2*(-(4*a*c - b^2)^3)^(1/2) + 10*a*b^5*c - 4*a*b^2*c*(-(4*a*c - b^2)^3)^(1/2))/(c^4*(4*a
*c - b^2)^3))^(1/3)*(b^4 + 2*a^2*c^2 - 4*a*b^2*c)*(b*(-(4*a*c - b^2)^3)^(1/2) + b^4 + 16*a^2*c^2 - 8*a*b^2*c))
/(8*c*(4*a*c - b^2)))*((3^(1/2)*1i)/2 + 1/2)*(-(b^4*(-(4*a*c - b^2)^3)^(1/2) - b^7 + 32*a^3*b*c^3 - 32*a^2*b^3
*c^2 + 2*a^2*c^2*(-(4*a*c - b^2)^3)^(1/2) + 10*a*b^5*c - 4*a*b^2*c*(-(4*a*c - b^2)^3)^(1/2))/(54*(64*a^3*c^7 -
b^6*c^4 + 12*a*b^4*c^5 - 48*a^2*b^2*c^6)))^(1/3)

________________________________________________________________________________________